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ABSTRACT 

Torsional  shafts  made  from  fibrous  composite  materials 
are  designed  for  minimum  weight.   Multiple  loading  conditions 
are  considered.   Design  variables  are  the  shaft  inner  diameter, 
the  ply  thicknesses,  and  the  ply  orientations.   Design  con- 
straints include  limits  on  ply  strain,  displacement,  frequency, 
and  Euler  and  shell  buckling.   The  design  task  is  solved  as  a 
numerical  optimization  problem.   Examples  are  presented  to 
demonstrate  the  method. 
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I.   INTRODUCTION 


Optimization  and  composite  materials  are  two  fields  that 
have  seen  considerable  developments  in  recent  years.   Advance- 
ments in  these  technologies  have  provided  the  design  engineer 
with  new  dimensions,  the  depths  of  which  are  yet  to  be  explored. 

The  major  material  characteristics  that  are  advantageously 
cultivated  in  composites  are  the  high  strength  and  stiffness 
to  weight  ratios.   The  weight  reduction  is  so  significant  that 
the  aircraft  and  aerospace  industries  have  concentrated  major 
research  efforts  in  this  field.   Investigations  in  improving 
the  resistance  of  composites  to  environmental  factors  such  as 
temperature,  corrosion,  and  wear  are  continuing  and  findings 
affirm  the  contention  that  the  future  of  composites  is  bright. 
However,  other  industries  are  slow  in  utilizing  these  materials 
in  basic  machinery  elements,  mainly  due  to  a  major  drawback; 
the  high  cost  of  manufacturing  and  processing.   It  is  in  this 
vein  that  the  importance  of  an  optimized  design  becomes  apparent. 
Advancement  in  fabrication  and  processing  of  composite  materials 
coupled  with  the  capability  to  optimize  designs  may  well  lead 
composites  into  a  more  competitive  market. 

The  application  of  numerical  optimization  techniques  to 
the  design  of  composite  shafts  is  demonstrated  here.   The  effects 
of  a  small  mass  imbalance  on  the  deflection  of  a  composite 
shaft  with  synchronous  whirl  are  included. 


Section  2  defines  the  scope  and  limitations  of  the  analysis 
and  optimization. 

Section  3  describes  the  analysis  of  composite  shafts, 
including : 

A.  Composite  properties 

B.  Shaft  section  properties 

C.  Composite  failure  criteria,  ply  strain  calculations 

D.  Displacement  and  vibration  analysis 

E.  Summary  of  failure  modes  considered 

Section  4  presents  the  optimization  results  [Refs.  1,  2]. 
Three  design  examples  are  shown. 

Section  5  contains  the  conclusions  drawn  from  this  study 
and  recommendations  for  future  investigations. 

Appendix  A  contains  a  description  of  the  analysis  portion 
of  the  FORTRAN  program  used  for  composite  drive  shaft  design. 
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II.   SCOPE  AND  LIMITATIONS 

The  analysis  presented  here  applies  to  the  shafts  made  of 
isotropic  or  composite  materials  with  pinned-pinned  end  con- 
ditions.  Basic  shaft  design  formulas  which  include  the  effects 
of  a  small  mass  imbalance  are  used.   The  shaft  is  assumed  to 
rotate  with  synchronous  whirl. 

The  optimization  is  constrained  to  subcritical  speeds  only. 
Supercritical  speeds  introduce  feasible  design  regions  that 
are  disjoint  from  the  primary  design  space.   In  this  investi- 
gation, these  secondary  areas  are  considered  infeasible. 

Since  the  optimization  tends  to  a  large  radius-small  thick- 
ness design,  radial  stresses  are  neglected.   However,  buckling 
as  a  thin  cylinder  due  to  torsion  or  compression  is  considered. 
The  equations  used  for  these  failure  criteria  do  not  include 
dynamic  effects.   A  constant  torque,  centrifugal  force,  and 
axial  force  are  considered  in  the  formulas  [Ref.  3]. 

The  analysis  incorporates   the  ability  to  design  a  shaft 
that  can  be  used  in  two  or  more  loading  conditions.    For 
example,  a  shaft  may  be  designed  to  transmit  50  HP  at  300  RPM 
with  an  axial  load  of  2000  lbs.  as  well  as  to  transmit  150  HP 
at  3000  RPM  with  a  1000  lbs.  axial  load. 
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III.   ANALYSIS 

The  problem  considered  here  is  a  shaft  required  to  transmit 
a  specified  horsepower  at  a  given  speed.   Basic  equations  are 
used  throughout  the  analysis. 

The  inch-pound-second  system  of  units  (IPS)  is  used  in  the 
computer  program  analysis.   However,  by  changing  the  equation 
for  torque  and  the  acceleration  due  to  gravity  to  reflect  the 
International  System  of  Units,  the  program  may  be  used  with 
consistent  SI  units. 

The  assumptions  used  in  the  analysis  are: 

1.  The  plys  are  made  of  orthotropic  material,  but  need 
not  all  be  the  same  material. 

2 .  S  true tur e : 

a.  Plane  stress  in  radial  direction 

b.  Plys  are  at  different  orientations  and  have 
different  thicknesses. 

c.  Cross-section  is  constant  over  entire  length 

3.  Loads:   (see  Figure  1) 

a.  Axial  force 

b.  Constant  torque 

c.  Centrifugal  force 

4.  Failure  modes: 

Typical  ply  orientation  and  composite  layup  are 
shown  in  Figure  2. 
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Figure  1 


Applied  Loads:   Axial  load  (F), 
Mass  imbalance  (e),  Torque  (T). 
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Figure  2.   Typical  ply  orientation  and  composite  layup 
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a.  Maximum  ply  strain  limits 

b.  Maximum  displacement 

c.  First  fundamental  frequency 

d.  Shell  buckling  (axial  and  torsional) 

e.  Euler  buckling  (contained  within  the  displace- 
ment constraint) 

A.   COMPOSITE  PROPERTIES 

Stress  is  a  measure  of  internal  force  within  a  body.   This 
together  with  elastic  constants  are  the  variables  for  the 
determination  of  stiffness  and  strength  of  a  material.   The 
mechanism  of  deformation  and  failure  are  interpreted  in  terms 
of  the  state  of  stress  and  strain.   When  discussing  stress  it 
is  usually  considered  to  be  the  average  stress  over  some  physi- 
cal dimension. 

In  Composites,  three  levels  of  average  stress  are  considered 

a.  Micromechanical  or  local  stress  is  that  calculation 
based  on  distinct,  continuous  phases  of  fiber, 
matrix  and,  in  some  cases,  the  interphase  and  voids. 

b.  Ply  stress  is  that  calculation  based  on  assumed 
homogeneity  within  each  ply  or  ply  group  where  the 
fiber  and  matrix  are  smeared  and  no  longer  recog- 
nized as  distinct  phases. 

c.  Laminate  stress  resultant  N  or  moment  M  is  based  on 
an  average  of  ply  stresses  across  the  thickness  of 

a  laminate.   The  individual  plies  are  smeared. 


The  coefficients,  or  material  constants,  of  the  stress- 
strain  relations  can  be  packaged  in  set  of  engineering  con- 
stants, compliance  components  or  modulus  components. 

a.  MODULUS  is  used  to  calculate  the  stress  from  strain 
This  is  the  basic  set  needed  for  the  stiffness  of 
multidirectional  laminates. 

b.  COMPLIANCE  is  used  to  calculate  the  strain  from 
stress.   This  is  the  set  needed  for  the  calculation 
of  engineering  constants. 

c.  ENGINEERING  CONSTANTS  are  the  carryover  from  the 
conventional  materials.   Designers  often  feel  more 
comfortable  working  with  the  familiar  engineering 
cons  tant s . 

d.  There  is  a  direct  relationship  between  the  modulus 
and  compliance.   One  is  the  inverse  of  the  other. 

The  stiffness  of  unidirectional  composites  is  governed  by 
the  same  stress-strain  relation  that  is  valid  for  isotropic 
materials.   However  the  number  of  independent  consitants  are 
four  for  planar  composites  as  compared  to  two  for  isotropic 
materials.   The  stiffness  of  isotropic  materials  can  be  repre- 
sented by  the  Young's  modulus  alone  because  the  Poisson's  ratio 
is  the  same  all  directions. 

For  composite  materials,  Poisson's  ratios  are  not  bounded 
and  can  have  very  significant  effect.   Young's  modulus  alone 
is  not  sufficient  to  describe  the  stiffness  of  composite 
materials.   It  is  not  sufficient  for  isotropic  materials  either 
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The  ply  stiffness  related  to  the  global  coordinate  system 
are  [Ref .  4] : 

/  L.  9  9  9   9 

Q,.,     =    m    Q         +nQ         +2mnQ         +4mnQ 
xll  xx  yy  xy  ss 

/  /  0      0  0      0 

Q22  =  n  Qxx  +  m  Qyy  +  2m  n  9xy  +  4m  n  Qss 

Q12  =  '^^n  +  m2n2qyy  +  ^^^^y   "  *»2n\8 

2  2         2  2  2  2   9. 
Q_  =  m  n  Q    +  m  n  Q    -  2m  n  Q    +  (m'-n  )Q 
x66        xxx         yy          xy  ss 

3  3  3   3  3   3 

Q,,  =  m  nQ    -  in  Q    +  (mn  -m  n)Q    +  2(mn  -m  n)Q„„ 
x16       xxx       xyy  xy  ss 

Q26  =  mn3Qxx  -  m3nQyy  +  (mn3-mn3)Qxy  +  2(m3n-mn3)Qgs 


where 


m  =  cos  9 


n  = 


xx 


sin  0 

E 
x 

1-1?  k> 

x  y 


•yy 


•xy 


l--v>  *> 
x  y 

V  E 

y  x 

x  y 


Q    =  E 
ss     s 


and  Q   ,0   ,  Q   ,0   :   On-axis  modulus  components  referred  to 
xx    yy    xy    ss  r 

the  local  coordinate  system 
Q.,:   Off-axis  modulus  components;  i,J  =  1,2,6 

1  U 


B.   SHAFT  SECTION  PROPERTIES 

Axial  stiffness,  bending  stiffness,  and  torsional  stiffness 
are  smeared  properties. 
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a.   Axial  stiffness  (AE) 


NPLY 
AE    =      SI      CQxl)kAk 
k=l 


(3B.1) 


b.   Bending  stiffness  (EI) 


NPLY 

ei  ■  s:  (Qn)kik 

k=l 


(3B.2) 


c.   Torsional  stiffness  (GJ) 

NPLY 


GJ    = 


k=l 


^llVk 


(3B.3) 


where 

h 

Jk 

Ak 

NPLY 


the  component  of  modulus 

the  ply  moment  of  inertia 

the  ply  polar  moment  of  inertia 

the  ply  cross-sectional  area 

the  number  of  composite  plys. 


i , j  =  1,2,6 
k  =  1,NPLY 
k  =  1,NPLY 
k  =  1,NPLY 


C.   COMPOSITE  FAILURE  CRITERIA 

For  the  determination  of  strength  of  any  material  it  is 
usual  practice  to  estimate  the  stress  at  the  time  and  location 
when  failure  occurs.   In  the  case  of  conventional  materials  it 
is  needed  only  to  determine  the  maximum  tensile,  compressive, 
or  shear  stress  or  the  principal  stresses  and  then  some  obser- 
vation can  be  made  about  the  failure  mechanism.   This  process 
is  relatively  straightforward  because  isotropic  materials  have 
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no  preferential  orientation  and  usually  one  strength  constant 
will  suffice.   The  isotropic  material  is  essentially  a  one- 
dimensional  or  one  constant  material.   The  Young's  modulus  for 
stiffness  will  suffice  because  Poisson's  ratio  is  the  same  in 
all  directions,  and  uniaxial  tensile  strength  will  also  suffice 
because  the  shear  strength  is  taken  to  be  about  50  to  60  percent 
of  the  tensile  strength. 

For  composite  materials,  however,  the  on-constant  approach 
for  stiffness  or  for  strength  is  no  longer  adequate.   Four 
elastic  constants  are  needed  for  the  stiffness  of  a  planer  ele- 
ment.  Six  constants  for  the  strength  of  unidirectional  composites 
are  needed.   Unidirectional  composites  have  direc t ionally   depen- 
dent    strengths.   The  longitudinal  strength  can  be  twenty 
times  that  of  the  transverse  and  shear  strengths;  so  for  any 
state  of  stress,  all  three  stress  components  must  be  examined 
before  a  judgement  on  the  cause  of  failure  can  be  made.   The 
specific  stress  component  that  is  responsible  for  the  failure 
may   not  be  easily  identified.   Probably  all  three  components 
are  responsible.   The  effect  of  combined  stresses  must  be  sys- 
tematically determined  and  can  be  regarded  as  a  way  of  life 
for  composites . 

For  composite  materials,  there  is  needed  a  failure  criter- 

ian  for  the  unidirectional  plies.   The  strength  of  a  laminated 

composite  will  be  based  on  the  strength  of  the  individual  plies 

within  a  laminate.   Successive  ply  failures  are  expected  as  the 

applied  load  to  a  laminate  increase.   The  first  ply  failure 

is  followed  by  other  ply  failures  until  the  last  ply  failure 
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which  would  be  the  ultimate  failure  of  the  laminate.   The  ply 
stress  and  ply  strain  calculations  are  intended  for  strength 
determination.   There  are  two  popular  approaches  for  failure 
criteria   of  unidirectional  composites.   They  are  each  based  on 
the  on-axis  stress  or  strain  as  the  basic  variable  [Ref.  4]. 

1 .  The  Maximum  Stress  Criteria 

^x  *  *x 

*"s  ^ 

Failure  occurs  when  one  of  the  qualities  is  met. 

2 .  The  Maximum  Strain  Criteria 

*•  x  —   x 

£s=  £"s 

Failure  occurs  when  one  of  the  equalities  is  met. 

The  maximum  stress  and  strain  criteria  are  not  the  same. 
Only  when  Poisson's  ratio  of  the  unidirectional  materials  is 
zero,  the  criteria  become  identical.   Conceptually  they  are 
similar.   Each  component  of  stress  or  strain  has  its  own  cri- 
terion and  is  not  affected  by  the  other  components.   There  is 
assumed  to  be  no  interaction. 

In  this  study  the  ply  strain  failure  criteria  was  used. 

D.   DISPLACEMENT  AND  VIBRATION  ANALYSIS 

In  order  to  do  the  displacement  and  vibration  analysis  the 

deflected  shape  of  beam  is  assumed  as: 
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0  =    A      sin  — p   (see  Figure  3) 
o      J? 


where  the  beam  is  assumed  to  be  simply  supported.   The  deflec' 
tion  is  maximum  at  x  =  1/2  [Ref .  3] 


/. 


max 


.  Hk 

K1e    +    K0 

384 

[    EI-Kl(|)*-F(|)2J 

(3D.1) 


where : 

Kx  =  ^Aw2 

K2  -/Ag 

^  =  mass  per  unit  volume 
A  =  cross  sectional  area 
»  =  shaft  speed  is  radians  per  second 
g  =  acceleration  due  to  gravity 
e  =  eccentricity  of  mass  with  respect  to  the  axis  of 

rotation 
■i  =  shaft  length 
EI  =  smeared  bending  stiffness 
F  =  axial  load 
Certain  characteristics  of  this  equation  are  of  particular 
interest.   If  there  is  no  rotation  or  axial  load,  the  formula 
reduces  to  the  classical  equation  for  maximum  beam  deflection 
under  its  own  weight. 


5K 


max    384EI 


(3D. 2) 


If  there  is  no  axial  load,  the  equation  becomes  unstable 

0   4 
(deflection  goes  infinity)  when  EI  =  £-,(#)    or,  solving  for 

the  critical  speed, 


&>   =  71" 

c 


EI 


r 


aJL4 


(3D. 3) 


which  defines  the  first  fundamental  frequency  of  a  simply 
supported  shaft. 

If  there  is  no  rotation,  instabilty  arises  when 


2 
77  EI 

J?2 


(3D. 4) 


recognizable  as  Euler's  column  buckling  criteria. 

With  both  rotation  and  axial  load,  the  critical  speed 
becomes 


»   =  77" 

c 


EI  -  F(£) 


Al4 


(3D. 5) 


To  avoid  the  instability  regions,  the  denominator  of 

equation  (3D.1)  must  be  given  a  lower  bound  greater  than  zero 

(point  A  in  Figure  4).   Additionally,  a  maximum  deflection  has 

to  be  imposed  to  prevent  computer  overflow  (point  3  in  Figure  4) 
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These  limits  have  no  relevance  with  the  deflection  constraint 
used  in  the  design  optimization  process  (point  C  in  Figure  4), 
but  are  used  only  to  prevent  numerical  ill-conditioning  on  the 
computer . 

E.   SUMMARY  OF  FAILURE  MODES  CONSIDERED 

Failure  criteria  serve  important  functions  in  the  design 
and  sizing  of  composite  laminates.   They  should  provide  a  con- 
venient framework  or  model  for  mathematical  operations.   The 
framework  should  remain  the  same  for  different  definitions  of 
failures,  such  as  the  ultimate  strength,  the  proportional  limit, 
yielding,  endurance  limit,  or  a  working  stress  based  on  design 
or  reliability  considerations.   Failures  in  composite  materials 
involve  many  modes,  such  as  fiber  failures,  matrix  failures, 
interfacial  failures,  delaminat ion ,  and  buckling.   Furthermore, 
the  various  modes  interact  and  can  occur  concurrently  or 
sequentially . 

NLC ( 6 CNPLY)+4)  modes  of  failure  are  considered  in  this 
investigation.   The  6(NPLY)(NLC)  give  the  maximum  strain  limits 
in  longitudinal  compressive  strain,  longitudinal  tensile  strain, 
transverse  compressive  strain,  transverse  tensile  strain,  nega- 
tive shear  strain,  and  positive  shear  strain.   Failure  is  assumed 
if  a  ply  strain  multiplied  by  a  safety  factor  is  greater  than 
the  corresponding   strain  limit. 

A  limit  on  the  deflection  is  imposed.   A  common  practice  is 
to  specify  a  maximum  deflection  in  inches  per   foot  of  shaft 


o  i 


«.  *^ 


*-  X 


Pinned-pinned  end  condition 
Figure  3.   Assumed  deflected  shape  of  the  shaft 


Deflection 


Denominator 


Figure  4.   Limits  on  the  deflection.   The  solid  curve 
represents  the  deflection  values  used  in 
the  analysis. 
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length.   For  example,  as  used  in  this  report, 


X  =  0.005  inch/foot  length 

v  max 

Since  the  optimization  may  lead  to  a  thin  cylinder  design, 
buckling  due  to  torsion  and  compression  are  considered.   The 
critical  torque  and  critical  compressive  stress  for  buckling  of 
thin  cylinders  are  [Refs.  5,  6]: 


775-    * 


T     ,     rrtf    e  ^rt5 

C         3(1-/)° 


_   **^Qii 


l/rr 


'"c    = 


Et 


tQ 


11 


r\/T(77)  ^r 


Experimental  data  indicate  that  actual  failure  of  the  cylinders 
usually  occurs  below  fifty  percent  of  the  values  calculated 
from  these  equations.   Therefore  a  large  factor  of  safety  is 
needed  here. 

Lastly,  failure  is  assumed  if  a  specified  percentage  of 
the  critical  speed  as  calculated  by  equation  (3D. 5)  is  less  than 
the  shaft  speed . 

The  factors  of  safety  used  in  this  study  are  listed  in 
the  following : 


FAILURE  MODE 

Strain 

Torsion  Buckling 

Shell  3uckimg 

Frequency 

Deflection 


SAFETY  FACTOR 
2 
5 
5 
5 
1 


IV.   OPTIMIZATION  RESULTS 

There  are  a  number  of  optimization  programs  available, 
each  one  using  different  technqiues  in  locating  the  desired 
optimum  design.   COPES/CONMIN  [Refs.  1,2]  is  a  versatile  pro- 
gram that  may  be  used  for  sensitivity  analysis  as  well  as  an 
optimization  tool.   Provided  with  a  user  supplied  analysis  pro- 
gram (Subroutine  ANALIZ)  where  the  objective  function,  con- 
straints and  other  relevant  parameters  are  calculated,  it 
determines  a  usable  and  feasible  sector  from  which  a  search 
direction  is  chosen.   This  choice  of  search  direction  is  made 
by  using  such  information  as  the  gradients  of  the  objective 
function  and  gradients  of  active  and  violated  constraints.   This 
iterative  process  of  minimizing/maximizing  the  objective  func- 
tion by  changing  the  design  variables  is  terminated  when  no 
further  improvement  can  be  made. 

In  this  study,  the  objective  function  to  be  minimized  is 
the  weight,  with  the  shaft  inner  diameter  and  the  ply  thicknesses 
and  orientations  as  design  variables.   The  analysis  subroutine 
incorporates  the  analysis  of  composite  shafts  described  in  the 
previous  sections.   Several  design  examples  are  shown  to  demon- 
strate the  effects  of  the  different  parameters  on  the  optimized 
design.   For  the  examples  considered  here,  the  loads  are  listed 
in  Table  1  together  with  the  optimum  weight  and  the  constraints 
are  listed  in  Table  2. 
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Table  1:   Optimization  Examples 


Axial 

Optimum 

No. 

HP 

RPM 

Load 
(lbs) 

Material 

Weight 
(lbs) 

1 

150 

300 

0 

Steel 

50.2 

2 

150 

300 

0 

G/E 

12.5 

3 

1000 

6000 

3000 

Steel 

26.5 

1500 

8000 

0 

4 

1000 

6000 

3000 

G/E 

10.8 

1500 

8000 

0 

5 

25000 

110 

40000 

Steel 

24822.0 

10000 

30 

5000 

6 

25000 

110 

40000 

G/E 

1872.5 

10000 

30 

5000 

Examples  No.  1  and  No.  2  are  designs  for  one  load  condition- 

Shaft  lengths  are  120  inches  here. 

Examples  No.  3  and  No.  4  are  designs  for  two  load  conditions. 

Shaft  lengths  are  60  Inches  here. 

Examples  No.  5  and  No.  6  are  designs  for  two  load  conditions. 

Shaft  lengths  are  240  inches  here. 

The  results  of  optimization  are  presented  in  Figures  5-16. 


Table  2:   Design  Constraints 


Number 


Cons  traint 


For  NLC  =  1 


1  through  6(NPLY) 

6(NPLY)+1 

6(NPLY)+2 

6(NPLY)+3 

6(NPLY)+4 


Maximum  Strain  Limits 

Buckling  due  to  Torque 

Buckling  due  to  Compression 

Speed 

Def lection 


For  NLC  =  2 

6(NPLY)+5  through  12(NPLY)+4 

12(NPLY)+5 

12(NPLY)+6 

12(NPLY)+7 

12(NPLY)+8 


Maximum  Strain  Limits 

Buckling  due  to  Torque 

Buckling  due  to  Compression 

Speed 

Def lee  t ion 


Etc .  for  NLC  >  2 
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COMPCSI'F     CPIVFSHAPT    OUTPUT 


NUMBER  OF 
NUMBER  PF 
NUMBER    C'.F 


PLYS 

VATF; 
LCAO 


ECCENTRIC ITY 


[AL     TYP^S       = 
CONDITIONS    = 


I 
1 
1 
O.IOOOOF  +  OO 


DIMENSIONS 

PLY  THlCKNFS<;  PEFCSNT  OIAMFTFr-  ThETA 

INSIDE   c.2o:oo=*oi 
1    C.IOOOOE+Ol       IOC. 00         0.40000c+0l         O.Q 


STIFFNFSS 

AF     =    Q.31071E+0S 

FI     =    0.'M87HcfO9 
GJ    »    0.27187F+09 


LOACS: 
L.C. 
I 


0.315L3E*C5 


0.0 


0.0 


H     P 

0,150 OOP *03 


:.-,oooo  r>0? 


LOAD   CGMO  ITICNI 
PLY  FFL 

1  0.0 


1 
3.F 


.      .  FPT 

u.w  100.00         0.0 

CRITICAL     SFrcr)  =    0.l627aE+03 

J     •  '    OEFLFCTICN    =    0.2*0^6-01 

■    0.31893P-t-0i 

=    0.1L310F+0* 


V.C  i   i   n.- 

MAXIMi^ 

WEIGHT 

VOLUME 


S.F. 
100.00 


FPLT  S.F. 

J.23132E-03  5.61 


Figure  5.   Initial  Design  for  Example  1 
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COMPHSITE     CPIVFSHAFT    OUTPUT 

NUMBER    OF     Pl.YS  =  1 

NUMrtC?    nF     MATERIAL     TYPC?        =  I 

NUM8F=    HP     LOAD    CCNJ  IT  ILNS    =  I 

ECCENTRICITY  =  0 .1O0OOE+0O 

DIMENSIONS 

PLY  T~UCKNCSS  cFRC'NT  DIAMETE?  THETA 

INSIDE    0.5o407F*01 
1    0.ai786E-0l        10U.G0         <3.53543E*01  0.0 

STIFFNESS 
AE    =    Q»4doS6F*08 
EI     =    0.2037CE+Orv 
GJ    =    0.  1<V2S9F+0S 

LOAOS: 

L    C  T  M  F  HP  <DV1 

I     *  0.31513E+05  0.0  0.0  oIl5DUOE*03  J. 30003 ~  ►  3  3 

LOAC    CHNO I TIGN  I 

PLY                 EPL  S.F.                 EPT                 S.F.                 =PLT                 S.F. 

I          0.0  100.00         0.0                           100.00         0.O4690E-03          2.  J I 

CRITICAL     SPEED  a    0.2971 7E+03 
MAXIMUM     DEFLECTION    =    0.73*0lE-02 

WEIGHT  =    0.50l9lc+0Z 

VOLUME  =    0.1 7793E+03 


Figure  6.   Optimum  Design  for  Example  1 
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COMPCS!T?     C51VFSHAFT    UJTPUT 


NUMBER    ^P     PLYS 

NUK3C«       c     v\T=OTAL     TY=>  =  :;  = 

NUMBER     If     L<~AO    CCNCITHNS  = 

ECCFNT^'CITY  = 


1 
5.100006+00 


DIMENSIONS 

PLY          THICKNCS3  PC»C6NT          0IAuET5o 
INS [OF    G.20000P*0l 

1  C.2SJ'-)CE«-00  2b. 00          0.,>  5  0  00  E+0  1 

2  C.25u0lif*O0  25. 00          J.  3000J  "+0  I 

3  0.2500CE+uG  25. GO         0.0500JF«Ol 
<*    C.25 OOGE*OJ  25. OU         0. -tOOQu  =*0  I 


"HETA 

0.0 

0.20000P*O2 

-0.2J J03=f 02 
0.90000=*0? 


STIFFN-SS 

AE    =    0.1  21  *"JB*09 

EI     =    0. L2  42tc+09 
GJ    =    O.J7225T+^d 


LOAOS: 
L.C. 
1 


0.31513E+05 


0.0 


J  .0 


H  .  0  .  j  o  v. 

0.  15000=f03  G.2000C5+03 


LOAC 

PLY 
1 
2 
3 

4 


cnr-n  iticn 

FPL 

O.^08ll=-03 
-<}.'♦  7b  136-03 

0.53l59r-C9 


CRITICAL     SPrcD 


MAX  I^'UN 

WEIGHT 
VOLUME 


DEFLATION    = 


I 

S.F. 

100. 0  0 

21.  00 

1  ti.OO 

10  0.00 
=    0 


ept  ; 

0.0 
-0.^081  lc- J3 
0.^761  3E-  )3 
-0.531596-09 
0662F>03 
0.  1^7^1c-0! 
0.633  3tF>02 
=    0.U31J6+O4 


100. 00 

^tj.13 

9.89 

100.00 


-°1  T  S. 

0. 105826-02 

0.^72735-03 

3.1134-3F-02 

■J.  16931E-02 


F. 
I  7.  1  ) 
18.92 
16-2'. 
1C.  a  J 


Figure    7 


Initial  Design  for  Example 


COMPGSIT-     CPIVCSH4FT    rjjTPUT 


NUM9EC  r'F 
NUNQeK  ?F 
NUMRCQ    nF 


ECCENT- [CITY 


PlYS 

|X4T=^[AL     TYP"        = 

LO'.n    CCM^ITIi'NS    = 


4 

L 

1 

0  .10  000 E*00 


DIMENSIONS 

PLY         TrilCKr.FSS  PEf»CEfJT          DIAMETER 
INS  10E    0.1  2bJ6c*Cl 

1  0.l^0lE-O7  0.00       o.i2eotc+oi 

2  0.1*0  725*00  50.00         0.16*215+01 

3  0.io372P*03  50. 3C          G.20035F+JI 
*    G.1V..G1E-U7  COO         0.2003ic*01 


TRET  A 

0.0 

0.32225-^32 

-.).32225c»-0? 

3.900005*02 


STIFFN— S 
AF     =    G.21642C*  38 
EI     =    0.764fclF* 07 
GJ     =    G.6  37d7F*07 


LOAOS: 
L.C. 
1 

T 
J.31513E 

*05 

LOAC 
PLY 

1 

2 

3       - 

4 

CrNO  IT  IGN 

"PL 
0.3 

0.L92C7E- 

-0.23435P- 

0.  L63UF- 

02 
J2 

0  6 

I 
S.F 

ICO. 

**  * 

3. 

1GJ. 

00 
<V6 
66 

0  0 

0.0  0.0 


EPT  S.F. 

0.0  100.00 

-0.192075-32  9.16 

0.23<.35E-J2  2.01 

-0.163115-18     100.00 
CRITICAL     SPF=0  =    0.1152^F*0J 

MAXIMO"     OEFLECTICN     =    O.iOC O^F-01 
WEIGHT  =    G.1253OF*02 

VOLUME  =    0.223  76E*03 


H     P 

0.15000F+03 


3  PM 

3.  30000- +03 


-oiT       S.F. 
3. 332065-02    5.c> 
3.10364E-G2   10.CZ 
J. 22*065-02    B.?' 
•3.5l950=-02    3.5- 


Figure  8.   Optimum  Design  for  Example  2 


COMPOSITE     CRIVFSHAFT    OUTPUT 


NUM3F=  HP 
NU'^F'  ">P 
NU^eFP    pP 


ECCENTP [CI TY 


PLYS  = 

C&T(=DIAL    TYPCS 
LOAD    CCNTIT  IONS    = 


O.IOOOOE  +  CO 


DIMENSIONS 

PLY  THICKNFS?  DPRCFNT  OIA*FTE<?  THE">"4 

INSIDE    C.10J30E*02 
1     C.lOOOOE  +  01        IOC. 30         0.12G00EK)2         0.0 


STIFFNCSS 
AF     =    Q.113S3E*  IC 
EI     =    0.17374E+11 
GJ    =    0.12U2E+11 


LOADS: 
L.C. 

I 

2 


T 

0.105046*05 

G.11817E+0b 


0.0 

0.0 


M 


0.300G0E-f04 
3.0 


0.  IQ0  00E+04 
0.  15000E+04 


90M 

r:.6  0000  =  +04 
0.-3000  J  r  +  04 


LOAC    CCNDI  TICN  1 

PLY  .    EFL  S.F.  EPT 

I  0.26333P-G5     ICC.OC         0.0 

CRITICAL     SPCED  =    0.22742E+04 

MAXIV'JM     DEFLECTION    =    0.10584E-01 


S.F. 

100.00 


C°LT  S« 

0.51323c-05 


F. 
100.  0  ; 


LOAC    CONDITION  2 

PLY                  EFL  S.F.                  EPT 

1          0.0  LOO. 00         0.0 

CRITICAL     SFPFQ  =     0.22743E  +  04 
MAXIMUM    DEFLECT!CN    =    0.2G022<=-01 

WEIGHT  =     0.58471  =  4-03 

VOLUVE  =    0.20735E+Q4 


S.F. 
100.00 


=PLT  S.F. 

0.58301E-35    100.00 


Figure  9.   Initial  Design  for  Example  3 
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COMPCSITF     CclVESHAFT    OUTPUT 


NUMB F 3    nF 
NUMRF9    rp 

NUVBfS    -= 


ECCEMPIC  !  TY 


FLYS 

MATERIAL     Tyo^       = 

ICAQ    CONDITIONS    = 


1 
I 

C  .10000 E+CO 


DIMENSIONS 

PLY         THICKNESS  PFP.CENT  OIAMFTE*  THETA 

INSIDE    G.99<?92E*01 
1    0.49525E-01        100.00  0.  10096 F  +  02  0.0 


STIFFNESS 
AE    =    0.51542E+08 
EI     =    0.65C59F>u9 
GJ    =    0.45541E+09 


LOADS: 
L.C. 
I 

2 


0.1G504E+05 
0.U817E+05 


j.0 
O.J 


0.30000S+OV 
0.3 


M.o. 

O.IOOOOE>04 
0.  150G0F*04 


PPM 

3.60003  r04 
C.S0000r+'J4 


LOAD    CONDITION  1 

PLY  =FL  S.F.  EPT 

1         0.5d2C5P-04       n.ld         0.0 
CRITICAL     SFEEO  =    0.206  73E+04 

MAXIMUM    DEFLECTION    =    C.13034E-01 


S.P. 

100.00 


=PLT  S.F. 

3.U646E-03        11.1 


LOAC    CONDITION  2 

PLY                 EPL  S.F.                  EPT 

1          0.3  1G0.O0         0.0 

CRITICAL     SPEED  =    0.20691F+04 
MAXIMUM    DEFLECTION    =    Q.25007?-01 

WEIGH1"  =    0.264546  +  02 

VOLUME  =    0.93807E+02 


S.F. 

100.00 


=PLT  S.F. 

3.13102E-33  9.9  2 


Figure  10.   Optimum  Design  for  Example  3 
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COMPCSITF     C?IVCSHAFT    CUTPUT 


r4UH8F  = 
NUMBFR 
NUMBtR 


C'F 


ECCENTP  re  I  7Y 


Fl  YS 

^Ax"  I-  I  AL    ~YP  =  S        = 

LHAfJ    CCNHIT  IPMS    = 


4 
I 
2 
O.IOOOOE+OO 


DIMENSIONS 

PLY          THlCK.\FSS  PFRCrNT           DIAMETER 
INS  !OE    C.60000E+01 

1  C.2500CF+00  25.00         3.e5000  -+0 I 

2  Q.25000F+00  2b. 00          0. 7GO00 r*G  I 

3  Q.2500GE+00  25. JO         0.7500U=*0l 

4  C.250QGE+00  2b. 00         0.  oGOOu  c-0  I 


THETA 

0.0 

0.20000=+02 

-0.20000c+02 

0.90CC0=+02 


STIFFNCSS 

AE    a    0.29855E+09 

EI    =    O.I  7  2  Ur*  10 

GJ    =    C.4DCeHc>G9 


LOAD 
L.C. 

L 

2 

5: 

T 
0.10504F*0^ 
0.11817F+05 

0 

0 

.0 
.0 

LHAC 
PLY 

1 

2 

3      - 

CONOI TICN 

FFL 

C.10C49C- 

0.35C6C=- 

-0.192C5F- 

0.29259=- 

0* 

04 
04 
10 

1 

S.F 

100. 
ICO. 
100. 
UO. 

00 

0  c 

0  0 
Oj 

- 

0.300  00E+0-* 
0.0 


CRITICAL     SF  =  cO 


=PT  S.c. 

0.0  100.00 

-0.7503  IE-04     110.00 

0.2925  4c-  J4     100.00 

J.10C49F-J4     100.00 


=     J.20liU=+04 


H      P 

oil  0000=*  04 
0.150  00=+04 


).  6  0000' +04 
C.S3Q0O~+O4 


_°LT  S.c. 

?.  7b7l6F-J4     10  0.  C: 
J.56005c-o4    100.0' 
0.  7  33  94E-04     100.0  j 
-0.  93  I  89  c- 04     IOO.Cj 


MAXIMUM    OEFLECTICN    =    0.1332SE-01 


LOAD    CCNO  ITICN             2 

* 

PLY                  EPL                  S.F.                   EPT                  S.P. 

EPIT                  S. 

F. 

1          0.0                           100.00          0.0                           100.00 

J. 35181F-04 

ioo. o: 

2          0.29^82E-04    100.00      -0.2948 2F- 04     100.00 

0. 7 02 72 E- 04 

100. OJ 

3       -0.3156tic-C4     IC0.00         0-3l5fcaF-04     ICO. 00 

0.75791 E-04 

I  00.0  '• 

4          0.32917C-10    100.00       -0.3291  7E-L0     100.00 

-J.  10484E-03 

100. 0. 

CRITICAL     SF^rn                =    j.^01 3*c+04 

MAXIMUM     OEFLECTICN    =    J.26695IS-01 

WEIGHT                                      =    0.7389UF+02 

VOLUME                                      =    Q.13195E+04 

Figure  11.   Initial  Design  for  Example  4 


COMPCSITP    CRTVESHAFT    nijTPUT 


NUMBER  r!c 
MJM8C-  ,f- 
NU^rtEP 


LP  ad 


ECCENTBICI TY 


CCN! 


TYP"        = 
IT  I  QMS    = 


4 
1 
2 

o.iooooe+oa 


DIMENSIONS 

PLY          ThICKNcSS  =>FRCcNJT          OTAMFTE' 
INSIDE    J.54431E+01 

1  0.1509<tE+U3  t'.6/  0.=T450E+01 

2  0.15O06E-OI  o.57         0.'>7762~+01 

3  0.15c36e-Gl  d.57         0.5807<f=-h01 

4  0.59~o7E-ll  (J. 00          Q.58074F*0l 


THETA 

0.0 
J.86400P+01 

-0.364005+01 
0.90000F+02 


ST!FFNCSS 

AE  =  O.^7 3  52  =+08 
EI  =  O.?o65«  =  +09 
GJ    =    O.labSS'+oa 


LOAOS: 
L.C. 

1 

2 


0.1C504E+G5 
0.118175+05 


0.0 
J.O 


0.30000E+04 

0.0 


H     P 

oIlGOOOE+04 
0.  l5000c+04 


3PM 

3.63000 =+04 

:.3  30oo  =  +O4 


ioac 

ply 

i 

2 
3 


CC'J  TITICN 

FFL 

0.44542=- 04 

0.2  8  5S0F-03 

-0.20CC3c-03 

0.5  I  49 lc- 09 


CITICAL     SP  =  cO 
MAXIMUM     "*EFLECTICN 


1 

S.F. 
I  CO. 00 

2=«.  19 

42.  o4 
ICO.  00 


=  PT  « 

0.0 
-0.2412  5F-03 
0.2445  7E-JJ 
0.44  542E-04 


=    0.2 06 71 =+04 
=    0.13J3wE-0l 


.F.  ^PLT                  S.F. 

130.00  3.16223F-02       11.34 

72.95  0.154<*3=-02        11.91 

19.25  3.15792F-02        U.G^ 

100.00  -U.16400E-02       11.22 


LOAC   CON 

0  I  T I  ON 

-> 

PLY 

FPL 

S.F. 

EPT 

S 

.F. 

I         0.0 

ICO. 

00 

0. 

0 

100. 

00 

2         0.2 

7254c-03 

31. 

44 

-0. 

2725-»E- 

03 

6t, 

58 

3       -0.274C2=-03 

31. 

23 

0. 

27402E- 

03 

17. 

19 

4          0.57928F-09 

100. 

Oo 

-0. 

5792 3E- 

o9 

100. 

00 

CRITICAL 

SF-=D 

= 

0. 

20  711' 

=  +04 

MAXIMUM 

OFFL^CTICN     = 

0. 

24942 

c-0l 

WE IGHT 

= 

0. 

108  1c 

=  +  02 

VOLUME 

= 

0. 

19314 

=+03 

=PLT  S.F. 

0.1S251F-02       10.03 
0.17522c-32       10.50 
3.176l7=-02       10.44 
-0.13450E-02  9.97 


Figure  12.   Optimum  Design  for  Example  4 
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COMPOSITE     CPIVESHAFT    O'JTPUT 


NUMB PR 

NUM8FR 
NUM8E' 
ECCENT? 


3C     FLYS 

OP     ViTcciAL     TYPrS 

D11    LPftD    CCNPITIONS 

ICITY 


I 
1 
2 
0.10000E  +  0O 


DIMENSIONS 

PLY         THICKNESS  PEPCFNT  OIAMETES  THETA 

INSIDE    C.5000CS+01 
I    0.80000E+01        100.00         0.2  10CC E+02         0.0 


STIFFNFSS 
AF    =    0.10771F+11 
EI     =    0.31 371P*12 
GJ    =    0.2!9fcCE+ 12 


LOADS: 
L.C. 

1 
2 


0.l<t3?*E-»-08 
0.21008E+08 


0.0 
0.0 


0.400COC*05 
0.500  00=  *o-v 


o!?5000F+05 

0.  10000c+05 


QPv( 

0.1  L000"O3 
0.  10000  c  +  02 


LOAO    CONDITION  I 

PLY  EPl  ".F.  EPT 

L  0.37136=-05     100.00         0.0 

CRITICAL     SPFcr)  =    o.  196366*03 

MAXIMUM     DEFLECTION    =    0. 131 79 E- 01 


S.F. 
100.00 


rPLT       S.F. 
0.68490E-03    l.^'J 


LOAC  CONO  ITION  2 

PLY       FFL  Z.F.                  CPT 

I          0.^642CF-06  ICO.OJ         0.0 

CRITICAL     SPCCT  =    J.196^3C>03 
MAXIMUM    DEFLECTION    =    0.12725E-01 

WEIGHT  =    0.22113E+05 

VOLUME  =    0.78414.E  +  05 


S.F. 

100.00 


EPLT  S.F. 

0.10045E-02  1.2* 


Figure  13.   Initial  Design  for  Example  5 
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COMPCS!TE     CPTVFSHAFT    OUTPUT 


NUMBFP    PF 
NUMBCP    OF 


ECCENTFICI TY 


PLYS 

VATE-IAL     "rYPc<i 

LOAD  CCNDITInNS 


1 
L 
2 
0  .10uOOE*00 


DIMENSIONS 

PLY  THICKNESS  PPRCpMT  DIAMETE3  THETA 

INSIDE    0. 120006*02 
1    0.63539E+01        100.00         0.24U8E+02  0.0 


STIFFNESS 
AE    =    0.120SLF+11 
El     *    G.57C52E+12 
GJ    =    0.39S26E+12 


LOADS: 
L.C. 

I 

2 


0.1<02<»F+C3 

C.21009E+C3 


0.0 
0.0 


0.400  00c+05 
0.50000?*0<t 


0.25000E+05 
0.  10000c+05 


-  PM 

O.I  1000" +03 

}.  30000' +02 


LOAD    CONDITION  I 

PLY  FPL  S.F.  EPT 

I         0.3  3033F-05    100.00         0.0 
CRITICAL     ^PE=0  =    0.2^99<iE+03 

MAXIMUM.  OEFLECTICN    =    0.81212E-02 


S  .F. 
100.00 


-PLT  S.F. 

0.^4328P-03  2.93 


LOAO    CQNO  I  TION  2 

PLY                  "L  s.f.                  cot 

I          0.M353F-06  ICO. 00          0.0 

CRITICAL     SPcEO  =    0.2500^  +  03 
MAXIMUM    OEFLcCTIGN    =    0.73531c-0? 

WEIGHT  =    0.2^322c+05 

VOLUME  =    0.38022E+05 


S.F. 
ICO. 00 


=  PLT  S.F. 

0.650l^c-03  2.0; 


Figure  14.   Optimum  Design  for  Example  5 
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COMPOSITE  CSIVESHAFT  OUTPUT 


NUM8FQ  nF 
NUMBER  Or 
NUMBER  PF 


ECCENTPICI TV 


PLYS 

f  4TERI  AL    TYPCS       = 
LOAO    CCNOIT  I  DNS    = 


L 

O.IOOOOE+OO 


OIMENSICNS 

PLY         THICKNESS  P=RCCNT          OIAMPTcR 
INSIOE    O.^OOOO^+Ol 

1  C.20C0CE+01  25.00         0.*0O00E+Oi 

2  Q.2000CE+01  25.00         0.13000F+02 

3  0.20000E+01  25.00         0 . 1 7000 F+02 

4  C.20000F+01  25.00         0.2 10CC c*02 


THETA 

0.0 

0.20000E*02 

-0.200  )0c+-02 

0.90CUOC«-02 


STIFFNESS 
AE    =    0.33749C*10 
El    =-  0.738S6F+U 
GJ    =    0.27562=*ll 


LOADS: 

L.C.                           T 

1  0.1<V324E*G8 

2  C.21006E+08 

M 
0.0 
0.0 

F                         HP 
0.^0000  =  +  05          o!;?5"!OOE*05 
0.50000  =  *OV          O.lOOOOE+05 

LOAO    CONDITION             1 
PLY                  EFL                  S.F 

1  0.  1032  3E-04    IGO. 

2  0.10948--D2          7. 

3  -0.141C6F-U2         f>. 
«•         0.17134—08     ICO. 

CRITICAL     SFc=0 

EOT 
00         0.0 
8  3       -0.l08-r5E- 
OH          0.L-+210E- 
00          0.1032  IE- 
0.22071E+03 

S.F.                  FOLT                  S.F. 

100.00         0.23387E-02           7.  a  7 

•02        16.23          O.?5811=-02          7.'3 

02          3.31          0.33r?06E-02          5.^i 

■J4     100.00       -0.5<V56>/P-02           3.~'7 

RPM 

0.1 100  >=♦•">  3 

0.3000;  :«-02 


MAXIMUM    OEFLECTICN    =    0.10<+22c-0l 


LOAO    CONDITICN            2 

PLY                  EPL                  S.F. 

EPT                  S.F. 

1          0.129C3F-05     ICO. 00 

0.0                           100.00 

2          0.1S935E-02          5.38 

-0.15922c-0>        11.05 

3       -0.20812F-02         *.12 

0.208Z5E-02          2.  ?o 

<+         0.25l29c-08     ICO. 00 

0.128^8E-05     100.00 

CRITICAL     SPFrD               =    0. 

221 02E+O3 

MAXIMUM    OEFLFCTICN    =    0. 

100  50E-01 

WEIGHT                                     =    0. 

<t39  l.^c+O^ 

VOLUME                                      =    0. 

78414E+05 

=PLT  S.F. 

0.3430l«=-02  5.36 

0.37946c-02  4.n? 

0.496*0E-02  3.7! 

-0.30035E-02  2.30 


Figure  15.   Initial  Design  for  ExamDle  6 


CCMPCSIT=     CRIVFSHAFT    OUTPUT 

NUMBER    r»P     FLYS  =  4 

NUMBER    OP     MAT^IAL     TYDCS  =  I 

NUMRPR    OP     LH4T    CC':n?TI°NS  =  2 

ECCENTRICITY  =  0.10003E+-00 

DIMENSIONS 

PLY         THICKNESS  PFRC=NT         OIAMETER  THETA 
INSIDE    0.10611*02 

1  0.123  5  I  F>00  3.b5         0.l085e=>02  0.0 

2  0.1542  7EOI  48.08         0.l3Q43eO2  0.42241  =  02 

3  0.154'  75*01  4«.oa         0.1702^=02  -0.42241FO2 

4  0.23842P-J6  0.00         0.17029EO2  0.90000EO2 

STIFFNESS 
AE    =    0.110165+10 
EI    =    0.271C9-+L1 
GJ    »    0^3778GE+H 


LOADS: 

L.C.                            T                                    M 

1  0.14324E+08         0.0 

2  C.21008EO8         0.0 

P                          HP 
0.40000EO5          ol25000E05 
0.50000FO'-»          Q.10000EO5 

LOAO    CCNO I TITN             1 

PLY                  PPL                  S.F.                  FPT 

1  O.36312E-04    ICO. 00          3.0 

2  0.13354E-02         6.42      -0.1299 1E- 

3  -0.15367P-02          5.40         0.16230F- 

4  0.  lOl3f.S-08    130.00         0.3631  IE- 
CPITICiL     SFEEO                =    0.l975dFO3 

S.P.                 CPLT                  S.F. 

100.00         3.20583F-02          8.9^ 

02       13.55         0.21305=03       84.3-1 

•02          2.90         0.34659P-03       53.09 

•04    100.30      -0.32232=02          5.70 

RPM 

3.11000*03 

0.3000  J  z+OP. 


MAXIMUM    DEFLECTION    =    0.13017E-01 


LOAO 

CONDITION 

2 

PLY 

FPL 

S.F. 

=PT 

S 

.F. 

I 

0. 

45389=- 

05 

IOC. 

00 

0. 

0 

100. 

00 

2 

0. 

19319=- 

02 

4. 

44 

-0. 

192  7  4F- 

•0? 

9. 

13 

3 

-0. 

23539P- 

0  2 

■3 

,64 

0. 

23584F- 

•0? 

2. 

00 

4 

0. 

14866=- 

08 

too! 

00 

0. 

453  75E- 

•05 

100. 

00 

CRITICAL     SF5=0 

s 

0. 

L9833F03 

MAXIMUM 

QEFL=C 

TION     - 

0. 

1243?=-01 

WEIGHT 

■= 

0. 

18725 

=  04 

VOLU 

*E 

= 

0. 

3343S 

£♦05 

EPLT  S.F. 

3.30139=02  6.09 

0.36831P-03       49. v6 
0.45984P-03       40.01 
-0.47347F-02  3.39 


Figure  16.   Optimum  Design  for  Example  6 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 


With  all  other  parameters  held  constant,  an  increase  in 

speed  drives  the  design  to  a  large  inside  diameter  while 
reducing  the  thickness  considerably.   This  trend  is  so  strong 
that  a  reasonable  lower  bound  on  the  thickness  must  be  given 
to  prevent  numerical  ill-conditioning. 

The  introduction  of  eccentricity  lets  the  optimization 
produce  larger  values  for  the  design  variables. 

The  effect  of  axial  load  is  reflected  on  the  deflection. 
If  very  small  deflections  are  allowed,  the  axial  load  effects 
are  reduced . 

The  versatility  of  the  program  should  not  be  lost  in  the 
emphasis  on  composite  shaft  applications.   The  objective  func- 
tion need  not  be  confined  to  the  weight.   If  the  outer  dimen- 
sion is  critical  to  the  design,  one  may  set  a  specific  value  on 
the  diameter  and  let  all  other  parameters  vary  during  optimization 

In  summary,  numerical  optimization  provides  an  efficient 
and  effective  way  of  creating  composite  shaft  designs. 

The  study  has  shown  the  feasibility  of  using  numerical 

optimization  techniques  in  the  design  of  composite  shafts  within 

the  limitations  imposed  in  the  analysis.   Further  studies  on 

the  same  design  field  may  be  pursued  by  eliminating  some  of  these 

limitations.   For  example,  the  unfeasible  region  defined  by 

the  speed  constraint  extends  only  to  some  discrete  distance 

and  a  supercritical  feasible  design  area  does  exist.   It  will 

hi 


be  interesting  to  develop  optimization  techniques  to  accomodate 
this  disjoint  design  field  and  to  investigate  the  peculiari- 
ties this  will  introduce  to  the  anlaysis. 

Another  investigation  may  be  concerned  with  shafts  with 
variable  thickness  and  diameter  along  its  length.   The  inclu- 
sion of  radial  stresses  in  the  calculations,  the  use  of  finite 
element  methods  in  three  dimensional  stress  analysis  of  the 
problem,  and  vibration  considerations  are  but  some  of  the  aspects 
of  research  that  may  be  useful  in  future  anlaysis  of  composite 
materials  . 


APPENDIX  A 


COMPUTER  CODE  USAGE 


The  FORTRAN  program  used  in  this  investigation  is  described 
here.   Subroutine  ANALIZ  provides  the  basic  analysis  used  in 
the  optimization.   It  reads  the  initial  design  descriptions  and 
calculates  the  values  of  the  objective  function,  constraints, 
and  all  other  parameters  necessary  to  solve  the   analysis 
problem.   COPES/CONMIN  updates  the  design  to  minimize/maximize 
the  objective  function,  iterating  until  no  further  improvement 
in  the  objective  function  is  possible  without  violating  one  of 
the  constraints.   This  process  is  shown  in  the  flow  diagram, 
Figure  A-l. 

The  global  catalog  given  in. Table  3  lists  the  locations, 
FORTRAN  name,  mathematical  symbol  and  description  of  the  para- 
meters used  in   the  optimization.   These  parameters  are  contained 
in  the  labeled  COMMON  block,  GLOBCM,  which  is  accessed  by 
COPES  for  optimization. 

There  are  seven  data  groups  to  describe  the  initial  design. 
The  first  data  group  contains  the  title.   One  card  is  used  here. 
The  second  data  group  contains  the  number  of  plies,  NPLY,  the 
number  of  different  materials,  NMT ,  and  the  number  of  load  con- 
ditions, NLC .   One  card  is  used,  here.   The  third  data  group 
contains  the  initial  inside  diameter,  shaft  length,  magnitude 
of  mass  imbalance,  fraction  of  shaft  radius  at  which  the  mass 
imbalance  acts,  and  the  acceleration  of  gravity.   One  card  is 

/  1 
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used  here.   The  fourth  data  group  describes  the  thicknesses 
and  the  orientations  of  plies  and  material  type  used  in  the 
analysis.   NPLY  cards  are  used  here.   The  fifth  data  group 
describes  the  ply  longitudinal  modulus,  EL,  ply  transverse 
modulus,  ET ,  ply  shear  modulus,  GLT,  ply  major  Poisson's  ratio, 
PRLT,  and  the  specific  weight,  RHO .   The  ply  transverse  Poisson's 
ratio,  PRTL,  is  calculated  internally.   NMT  cards  are  used  here. 
The  sixth  datagroup  specifies  the  ply  longitudinal  compressive 
strain  limit,  EPLC,  ply  longitudinal  tensile  strain  limit, 
EPLT,  ply  transverse  compressive  strain  limit,  EPTC ,  ply  trans- 
verse tensile  strainlimit,  EPTT,  and  ply  maximum  shear  strain 
limit,  GMLT  .   The  compressive  strain  limits  are  negative  num- 
bers, the  tensile  and  the  maximum  shear  strain  limits  are 
positive  numbers.   NMT  cards  are  used  here.   The  seventh  data 
group  describes  the  torque,  moment,  axial  load,  horsepower, 
and  shaft  speed  in  RPM  for  each  loading  condition.   NLC  cards 
are  used  here. 

The  horsepower,  torque,  and  shaft  speed  in  RPM  are  inter- 
dependent values.   If  two  of  them  are  known,  the  third  one  is 
calculated  internally.   The  following  table  summarizes  the 
required  information  for  Subroutine  ANALIZ. 
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Table  3:   Input  Data 


Data 
Group 

1 
2 
3 
4 
5 
6 
7 


Information 

Title 

NPLY,  NMT,  NLC 

DIAMI,  AL,  EO,  EC,  GRAV 

TN(I) ,  THND(I) ,  MTYP(I) 

EL,  ET,  GLT,  PRLT,  RHO 


Format 

2  0A4 
8110 
8F10.2 
2F10. 2,110 
8F10.2 


EPLC,  EPLT,  EPTC,  EPTT,  GMLT   8F10.2 
T,  M,  F,  HP,  RPM  8F10.2 
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SUBROUTINE  ANALIZ  (ICALC) 
Flow  Diagram 

ICALC  =  1 

Call  Subroutine  SHFT06 

Read  Initial  Dimensions 
Read  Material  Properties 
Read  Design  Conditions 
Print  all  Input 

ICALC  =  2 

Calculate  Geometric  Parameters 
Call  Subroutine  SHFT07 

Calculate  Objective  Function 
Call  Subroutine  SHFT03 

Calculate  Shaft  Section  Properties 
Call  Subroutine  SHFT09 

Calculate  All  Constraints 
Call  Subroutine  SHFT10 

Calculate  Critical  Frequency 
Call  Subroutine  SHFT11 

Calculate  Maximum  Deflection 

ICALC  =  3 

Call  Subroutine  SHFT08 

Print  Analysis  Results 
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OPTIMIZATION  FLOW  DIAGRAM 


INPUT 

' 

t 

COPES 

CONMIN 

< 

i       i 

' 

SUBROUTINE 
ANAL I Z 

\ 

> 

OUTPUT 

Figure  A-l 
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Table  4:   Global  Catalog  of  Parameters 

Global            FORTRAN  MATH 

Location           Name  Symbol        Definition 

1  VOL  V            Volume 

2  DIAMI  d           Inner  Diameter 

3  WGHT  w           Weight 

4-10  Dummy  Storage 

11-30            TN(20)  t           Thickness 

31-50            THND(20)  9           Orientation 

(Degree) 

51-551           G(500)  G           Constraints 
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